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ABSTRACT 

We present results of XMM-Newton Guaranteed Time observations of the RS CVn binary Coronae Bore- 
alis. The spectra obtained with the Reflection Grating Spectrometers and the European Photon Imaging Camera 
MOS2 were simultaneously fitted with collisional ionization equilibrium plasma models to determine coronal 
abundances of various elements. Contrary to the solar first ionization potential (FIP) effect in which elements 
with a low FIP are overabundant in the corona compared to the solar photosphere, and contrary to the "inverse" 
FIP effect observed in several active RS CVn binaries, coronal abundance ratios in cr^ CrB show a complex pat- 
tern as supp orted by similar findings in the Chandra HETGS analysis of cr^ CrB with a different methodology 
iOsten et al.li2003.) . Low-FIP elements (< 10 eV) have their abundance ratios relative to Fe consistent with the 
solar photospheric ratios, whereas high-FIP elements have their abundance ratios increase with increasing FIP. 
We find that the coronal Fe abundance is consistent with the stellar photospheric value, indicating that there is 
no metal depletion in cr^ CrB. However, we obtain a higher Fe absolute abundance than in Osten et al. (20 031). 
Excep t for Ar and S, our absolute abundances are about 1.5 times larger than those reported by Oste n et alJ 
However, a comparison of their model with our XMM-Newton data (and vice versa) shows that both 
models work adequately in general. We find, therefore, no preference for one methodology over the other to 
derive coronal abundances. Despite the systematic discrepancy in absolute abundances, our abundance ratios 
are very close to those obtained by Osten et al. (2003). Finally, we confirm the measurement of a low den- 
sity in O VII (< 4 X 10'" cm"^), but could not confirm the higher densities measured in spectral lines formed at 
higher temperatures derived by other studies of cr^ CrB due to the lower spectral resolution of the XMM-Newton 
grating spectrometers. 

Subject headings: stars: activity — stars: coronae — stars: individual {cp- CrB) — stars: flare — stars: late-type 
— X-rays: stars 



1. INTRODUCTION 

Studies of stellar coronae with X-ray missions such as 
XMM-Newton and Chandra have shown peculiar abundance 
patterns when compared to the the coronal composition of 
the Sun. In the latter, elements with a first ionization po- 
tential (FIP) below about 10 eV appear, on average, over- 
abundant compared to the high-FIP elements by a factor 
of a few (e.g., Mever 1985). In addition, high-FIP ele- 
ments are_of_solar photosphe ric comp osition (e.g., Feldma^ 
ll992HLaming. Drake. & Widingll995h . This abundance pat- 
tern is known as the "solar FIP effect". Previous studies 
of stellar coronae with lower-resolution detectors than avail- 
able onboard with XMM-Newton and Chandra and in the ex- 
treme ultraviolet rang e have shown a general depletion of the 
metal abundance (e.g.,'Antunes et al. fgg^.'White et al Jl994t 
ISchmittet al. 1996), and evidence for a solar-like FIP effect 
in some stars (Drake et al. 1997; La rning & Drakell999l) or its 
absence in others ( Drake et al. 1995^ 

XMM-Newton and Chandra have reached unprecedented 
spectral resolution and sensitivity in the X-ray regime thanks 
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to their grating spectrometers. This l ed to the discovery of 
an "inverse" FIP effect in HR 1099 d Brinkman etani200U 
Drake et al. 2001), in which low-FIP elements are under- 
abundant compared to the high-FIP elements. Other studies 
have strengthened the evidence for the presence of an inverse 
FIP eff ect in magnetically ac ti ve stars (e.g., ^ udard et afl 
'200 la* "Giidel et al.' '2001a'b' 'Huenemoerder et al. 20oJ 
2003; Audard etal. 2003; van den Besselaar et al. 20031 
IRaassen et al. 2003; Ost en et al.l i2003: Sanz-Forcad a et al] 
l200 3a). but have also shown that intermediately active stars 
show no strong correlation with the FIP (Audard et al. 2001^, 
and less active stars show a solar-like FIP effect ( Giidel et all 
|2P02; Telleschi et al. 2005). A transition from the inverse to 
the solar-like FIP effect thus seems to occur with decreas- 
ing magnetic activity (Giidel et al. 2002; Audard et al. 2003). 
However, the inactive F star Procyon does not show any cor- 
relation between its coronal abundances and the elemental 
FIP (Drake et al. 1997; Raassen et al. 2002), although its cool 
corona (1-2 MK) should display a strong solar-like FIP ef- 
fect in the transitional picture. Sanz-Forcada et al. (2004) ar- 
gued that the comparison of coronal abundances to solar pho- 
tospheric abundances instead of stellar photospheric abun- 
dances (which are often unknown or uncertain in the litera- 
ture) could be misleading as it could falsely mimic inverse or 
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solar-like FIP effects in stars. On the other hand, the pho- 
tospheric composition of the solar analogs in the study of 
iTelleschi et a is known to be similar to the Sun's; 

therefore, it appears that coronal abundance anomalies occur, 
at least in their sample. 

Several theoretical studies have p roduced m odels to explain 
the solar FIP effect (see reviews bvlHenouxll 995. 1998, also 
lArge et alJl998t lSchwadron et al.'1999^.'Mc Kenzie.200a) . In 
the stellar context, Giidel et al. (2002) suggested that the non- 
thermal electrons observed in the radio could explain the in- 
verse FIP effect in active stars, and possibly the FIP effect in 
inactive stars. Recently, Laming (2004) proposed a model that 
unifies the FIP and inverse FIP effects by exploring the effects 
on the uppe r chromo s p heric plasma of the wave ponderomo- 
tive forces. Lamina (2004) suggested that the observed so- 
lar FIP effect could be turned into an inverse FIP effect af- 
ter fine tuning model parameters. Further theoretical works 
are needed to address the wide range of FlP-dependent coro- 
nal abundances in the Sun and stars. In the mean time, we 
need to increase the sample of stars and to study their coronal 
abundances to understand the mechanism at the origin of the 
element fractionation in the Sun, and in stars. In this paper, 
we present the XMM-Newton study of the RS CVn binary cr^ 
Coronae Borealis which is composed of two solar analogs as 
well, albeit spun up by tidal forces. 

2. THE RS CVn BINARY a- CORONAE BOREALIS 

CrB (HD 146361; TZ CrB) is an active RS CVn bi- 
nary system at a distance of 21.7 pc ([Perrvmanet al. 1997) 
whose components have spectral types of F9 V and GO V 
((Strassmeier & Rice 2003). CrB is a component of a vi- 
sual binary with an orbital period of 1000 years, with the 
other component (cr' CrB) at a separation of 6.6". No eclipse 
takes place in CrB, due to the low inclination angle of 28° 
(|Bakos"1984), and the orbit has a period of P = 1.14 days 
(e.g., Drake et al. 1989). Strassmeier & Rice (2003) reported 
a 0.017 d difference between the rotation and orbital periods 
and interpreted it as differential rotation of the surface. The 
masses and radii of the com ponents are close to solar, i.e ., 
M ~ 1 . IMq and 7? ~ 1 .2Rq (|Barden'1985';"Drake et al."1989'). 
A detailed hst of the physic al parameters of^ cr^ CrB can 
be fou nd in lStrassmeier etal] l fT993h and lStrassmeier & Ricd 

(Hool. 

The binary CrB has been ob s erved extensively from 
radio to X-rays (see "Oste n et al] 12000*. and references 
therein). In the X-ray and EUV regimes, several studies 
have invest igated the properties o f the coronae of CrB 
dAgrawal et al 1981. 1 98 5. A9M Ivan den Oord 67^(1981 
lPasauinietalTl9 89: Ster n et alJl992HOsten et alJ2000h . Re- 
centlv. lOsteiiet'alj (i2003i) presented the Chandra spectrum of 
(7^ CrB obtained with the High-Energy Transmission Grat- 
ing Spectrometer (HETGS) as part of a coordinated observing 
campaign in the radio, EUV, and X-ray regimes. They found a 
bimodal emission measure distribution with peaks at 6-8 MK 
and 30 MK. The Fe abundance was found to be depleted 
as fre quently seen in bright RS CVn binaries ( Audard et ^ 
120031) . whereas Ne/Fe and Ar/Fe abundance ratios were above 
solar ratios. Although these ratios are indicative of an inverse 
FIP effect in CrB, other elements did not show a specific 
correlation with the FIP. In particular, the low-FIP (6 eV) ele- 
ment Al was overabundant with respect to Fe. With our anal- 
ysis of the XMM-Newton spectra, we aim to compare derived 
abundances with those of Osten et al- (2003). We will discuss 
the coronal abundances of CrB in the context of the FIP 



and inverse FIP effects in the Sun and stars. 

3. OBSERVATIONS AND DATA REDUCTION 

The o bservations of C rB were performed by XMM- 
Newton dJansen et alJ 120011) and were spread over three 
epochs: 2001 August 29 (hereafter "101"), 2001 August 31 
("201"), and 2002 February 22/23 ("301"). The first obser- 
vation lasted for approximately 7 ks and was interrupted by 
solar flare activity, therefore we will not discuss this data set 
further (except as a light curve in Fig.^. Solar activity was 
again high for the two remaining epochs; nevertheless, we 
obtained spectra with sufficient quality. We provide a log of 
these observations in Tabled 

All data were reduced with the XMM-Newton Science Anal- 
ysis System (SAS) 5.4. 1 and the calibration files of June 2003 
using standard procedures (e.g., Audard et al. 2003). The Eu- 
ropean Photon Imaging Camera (EPIC; Turner etal. 200ll) 
MOS2 data (in SMALL WINDOW mode) were shghtly 
piled-up, and therefore we excised the center o f the point- 
spread function (PSF) to avoid this problem (e.g.. Lum b et alJ 
2000). A background template was obtained from a region on 
an outer EPIC MOS2 CCD. Its exposure was corrected for vi- 
gnetting since we assumed here that the dominant X-ray back- 
ground was due to the cosmic X-ray background (i.e., the in- 
strumental and particle backgrounds, which are not vignetted, 
were neghgible). The EPIC MOSl data (in TIMING mode) 
were not used due to the small window size which prevented 
us from defining an appropriate background in this mode. The 
EPIC pn data were taken in TIMING mode as well and will 
also not be discussed here. 

Despite the high solar activity at the end of each observation 
(FigO, w e used the comp lete Reflection Grating Spectrom- 
eter (RGS: Iden Herder"et al...2001.) data, since a comparison 
with the "cleaned" data (selecting periods with low solar ac- 
tivity) showed that they were similar However, we preferred 
to use the cleaned EPIC spectra since they showed no sig- 
nificant contamination at high energies compared to the full- 
exposure spectra. We show the selected good time intervals 
(GTIs) in Fig.[l] Thanks to the low count rate variability of 
(7^ CrB during the XMM-Newton observations, the different 
treatments of GTIs between the EPIC and the RGS data had 
negligible impact. 

In our EPIC images, there was no evidence for X-rays 
from the visual binary component ct' CrB. Its X-ray emis- 
sion could in principle be separated (6.6") if it was bright 
enough (indeed, both Castor AB low-mass components were 
detec ted wifli XM M-Newton EPIC MOS for a separation of 
3.9": lGudel et al.^2001b:) . It suggests that the X-ray luminos- 
ity of (7^ CrB is much fainter than that of CrB. The Chan- 
dra zeroth order image confirms the faintness the component 
(R. Osten 2003, priv. comm.). Therefore, no significant con- 
tamination by (t' CrB is expected in the EPIC and RGS data 
of 0-2 CrB. 

4. DATA ANALYSIS 

The reduced RGS and EPIC MOS2 spectra were fitted si- 
multaneously using collisional ionization equilibrium (CIE) 
plasma models with variable elemental abundances and tem- 
peratures. A free constant multiplicative model was allowed 
to deal with cross-calibration normalization effects between 
the R GS and the EPIC spectra. W e used X SPEC (Arnaul 
119961) vl 1.2.0 softwai-e with APEC (Smith et al. 2001) code 
vl . 3.0 that contains line and c ontinuum emissivities. We used 
the|G revesse & Sauvall (Il998h set of solar photospheric abun- 
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dances. Such a set was also used bv lOsten etal] ( l200l thus 
allowing a direct comparison between their and our results. 

We used a multi-T approach with free temperatures and 
emission measures (EMs), and variable abundances (C, N, 
O, Ne, Mg, Si, S, Ar, Ca, Fe). We added a photoelec- 
tric absorption com ponent using Wisconsin cross-sections 
of Morrison & McC ammonI j 19831) but le ft the value of the 
column density fixed at 2 x lO'^ cm"^ (l Osten et aLll2000h 
since t he XMM-Newton da t a are n ot sensitive to such a low 
value. Sanz-Forcada et al. (2003b) suggest a value of 2.5 x 
10"^ cm"^, but the difference is negligible. We fixed the abun- 
dance of Ni to be equal to the abundance of Fe. 

Furthermore, we used a similar approach as lAudard et al.l 
(|2003). Several wavelength ranges in the spectra were dis- 
carded since they include low-Z L-shell transitions which 
are insufficiently described in current atomic databases 
(Q^udard et al. 2001b, 2003; Lepson et al. 2003). With this 
procedure, we have derived the element abundances from the 
more reliable emission lines of low-Z K-shell ions. The full 
list of exc luded w avelength ranges is given in Table|2] 

Figures l2al2bl show the EPIC MOS2 and RGS spectra over- 
laid with a best fit model using a 4-temperature CIE model. 
The EPIC M0S2 data at wavelengths shorter than 18 A were 
used to model the high-energy tail of the X-ray spectrum and 
the emission lines inaccessible to the RGS; significant over- 
lap with the RGS (8 - 18 A) was kept to make sure the high-T 
component is well linked to the lower- T components which 
mainly produce emission lines in the RGS band. The hot 
Fe Ka complex at 6.7 keV (weak in the EPIC MOS2 data 
but clearly detected in the EPIC pn data not shown here) in- 
dicates the presence of hot plasma whereas the O VII triplet 
reveals cool plasma. We list the best-fit parameters of A-T 
CIE models for the 201 and 301 observations in Table |3l We 
also list the abundances reported by Osten et al. ( 2003) for the 
Chandra observation of cr^ CrB. 

5. DISCUSSION 
5.1. Coronal Abundances 

We derived absolute abundances (i.e., abundances rela- 
tive to H ) with respect t o the solar photospheric abundances 
given bv lGrevesse & Sauv al ( 1998). However, abundance ra- 
tios, e.g., relative to Fe, have been found to be more robust 
(e.g., [Audard et al. 2004; Telleschi et al. 2005). Therefore 
we give abundance ratios to Fe in Table |3l Figure |3] shows 
the abundance ratios as a function of the FIP using the solar 
photospheric set of Grevesse & Sauval ( 1998). We plot the 
abundance ratios obtained for the 301 observation shifted by 
+0.2 eV for clarity. The figure shows two regimes: below 
10 eV, low-FIP elements have about the same abundances and 
are close to the solar photospheric values. Above 10 eV, the 
abundance ratios increase with increas ing FIP, as observed 
in sever al other RS CVn bi naries (Audard et alJl2003) . Fig- 
ure 1 1 o f lOsten et alJ ( l200l shows a similar pattern, although 
these authors preferred to interpret it as an absence of FIP- 
dependent pattern. 

Figure |3 shows a comparison of coronal abundances ob- 
tained in the best fits to the XMM-Newton 201 and 301 obser- 
vations, a nd extends it to a c omparison of the abundances ob- 
tained bv* Osten et alJ ( l2003h . It shows that abundance ratios 
are very similar for the three sets. On the other hand, the abso- 
lute abundances are systematically higher in our fits than those 
bv Osten et al. (2003.) by a factor of about 1.5. Although two 
different observations (e.g., at two different epochs) could, in 



principle, yield different abundances because the dominant X- 
ray emitting regions are not the same and the observed, av- 
erage, coronal composition could be different, it is far more 
probable that different analysis techniques could lead to dif- 
ferent best-fit abundances (e.g., Audard et al. 2004). Further- 
more, the spectral inversion problem being an ill-posed prob- 
lem, small statistical va riations can lead to discr epant results 
via the DEM inversion dCraig & Brownll976allJ) . It is, there- 
fore, difficult to estimate which set of abundances is preferred. 
On the other hand, it is worthwhile to stress the robustness of 
abundance ratios which show (whatever the technique and the 
data set) that a complex FlP-dependent abundance pattern is 
present in cr^ CrB's corona. 

The use of a solar photospheric set of abundances to com- 
pare with stellar coronal abundances is naturally problematic, 
since stars can have a photospheric composition at variance 
from the Sun's. Magnetically active stars are difficult ob- 
jects to derive accurate surface abundances in view of their 
rapid rotation and surface spots. Strassmeier &Rice (200t 
have, however, obtained estimates of the photospheric abun- 
dances of the low-FIP elements Fe and Ca, Ape = 7.30 and 
Aca = 6.07 (with the usual Ah = 12 notation; they estimate an 
uncertainty of ^ 0.1 dex ). Comparing with the solar values 
from Grevesse & Sauvall (|l998) (Ape = 7.50 and Aca = 6.36), 
(7^ CrB's photosphere appears depleted with [Fe/H] =-0.20 
and [Ca/H] =-0.29. Our spectral fits obtain [Fe/H] - -0.15 
which is consistent with th e stellar photospheric value. In 
contrast, 'O sten et alJ ( l200l obtain a slightly depleted coro- 
nal Fe abundance. It seems safe to conclude that the coronal 
Fe abundance in cr^ CrB is close to the photospheric value, in 
contrast to the solar case where the coronal Fe abundance is 
enhanced by a factor of a few. This situation is also different 
from that in several magnetically active stars that show a sim- 
ilar U-shape abundance pattern but a net Fe depletion in their 
corona (e.g., 'Gudel e t alJl200Ti |Sanz-Forcada et al." '20034 
l^lleschi et al. 2005). It is difficult to discuss the case of Ca 
since we derived large uncertainties and Os ten et alJ d20()3l) 
did not obtain a Ca abundance. Nevertheless, a scenario ap- 
pears to emerge in which low-FIP elements are depleted in 
very active stars that show an inverse FIP effect. The abun- 
dances of the lowest-FIP elements in intermediately active 
stars rise up to be quasi-photospheric while the intermediate- 
FIP elements still have their abundances depleted, thus the 
abundance pattern shows a U shape. Finally, in the less ac- 
tive stars, low-FIP elements have their abundances further in- 
creased to become superphotospheric and intermediate-FIP 
elements have their abundances quasi-photospheric. Such 
stars show an abundance pattern similar to the solar FIP ef- 
fect. 

5.2. A Comparison with the Osten et al. Model 

Our methodological approach to the XMM-Newton data 
(i.e., multi-T components) differs significantly from the ap- 
proach by Osten et al. (2003) (i.e., separate treatment of 
continuum and lines, reconstruction of a differential emis- 
sion measure [DEM] distribution based on extracted Fe line 
fluxes). We focus this section on comparing their model with 
ours. Specifically, we test whether their Chandra model can 
provide an adequate fit to our XMM-Newton dat a, and vice 
versa. Some caveats should be mentioned: i) Oste n et all 
(2003) use the APEC 1.10 version whereas we use the more 
recent APEC 1.3 version^; ii) CrB's X-ray luminosity 

^ http://cxc.harvai'd.edu/atomdb/index.htmI for a detailed comparison. 
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was 1.39 times higher during the XMM-Newton observation 
(Lx = 3.8 X 10-''' erg s"') than during the Chandra observation 
(Lx = 2.7 X 10^"' erg s~') . In our comparison, we will take 
this correction factor into account. However, it remains un- 
clear whether this factor should be applied to all temperature 
components or only part of them, iii) It is possible that real 
changes in the X-ray spectrum (e.g., due to EM or abundance 
variations) occurred between the two observations. We will 
assume in this study that this was not the case. 

We have used Osten et al.'s quiescent DEM (realization #2) 
together with the quiescent abundances reported in their Table 
5 (we set the C and Ca abundances equal to the Fe abundance, 
and we set Ni=0. 1 times solar photospheric, as described by 
[Osten et al. 2003). However, after iterations with Dr. Osten, 
we discovered that their quiescent DEM (as shown in their 
Fig. 10a) is accurate from log T{K) = 6.0 to 7.4 for the contin- 
uum and that the last point at log T{K) = 7.5 was used only for 
the emission lines. However, this temperature bin essentially 
contributed a few tens of percent to the flux in the Fe Ka com- 
plex. This oversight plays no significant role for our compar- 
ison exercise with the XMM-Newton data. Thus, in practice, 
we constructed in XSPEC a multi-T model with 15 compo- 
nents at temperatures at log T(K) = 6.0 to 7.4 with intervals of 
A log r = 0.1 dex. We used Osten et al.'s DEM values, multi- 
plied them by the original DEM bin width (A log T = 0.1 dex) 
to obtain the total EM for each temperature component, and 
further divided by 2 since Osten et al. (2003) accounted for 
the fact that half the radiation is absorbed by the star. 

We have then convolved the model of Oste n et alJ (120031) in 
XSPEC^ (using our RGS RMF file); we have also convolved 
our XMM-Newton EPIC M0S2 (301) model through the 
ACIS-S/HET GS first -order RMF/ARF provided by Dr Os- 
ten^. Figures I5al5cl compare the XMM-Newton RGS data 
overl aid wi th Osten et al.'s Chandra model, whereas Fig- 
ures I6al6cl show the Chandra HETGS data (provided by 
Dr Osten) overlaid with the XMM-Newton model. 

We find surprisingly good agreement, despite different ab- 
solute abundances. There are, however, some discrepancies. 
In Fig|5b] the Fe XVII A15 A line is overestimated compared 
to the XMM-Newton data. But this is no different from the 
analysis of Osten et al. (2003) which overestimated the Chan- 
dra Fe XVII line flux as well ( see their Fig. 10a, center panel). 
In addition, the lOsten et aP model overestimates the RGS 
continuum level by about 40% at wavelengths longer than 
25 A (Fig. l5ct . A few possible explanations can be put for- 
ward: i) the Chandra HETGS wavelength range ends around 
25 A, therefore the lOsten et all model had no possibility to 
accurately take into account the flux beyond this instrumen- 
tal limit. ISako et al. (2003) reported a similar situation in 
which the model of ILee et al.l (2001 ). based on HETGS data, 
could reproduce the general properties of the RGS spectrum 
of MCG -6-30-15 below 23 A, but it overestimated the flux 
for A > 25 A; ii) calibration problems at long wavelengths 
are possible; however, the 40% overestimate is larger than 

' IQsten et alJ |2003D used the composite trapezoidal rule to integrate the 
line and continuum flux at a specific wavelength. The composite trapezoidal 
rule helps to calculate the integral of a function f(x) over an interval [a, b] 
subdivided into N subintervals [.v,, of widths h = (b — a)/N (thus Xj = 
a + i X h for ; = 0, 1, . . . ,N). The integral is then equal to A X (f(a)+f{b))/2+ 
h X Y,f=i^ fM- Thus, we used half the DEM values at log T{K) = 6.0 and 
7.4 in XSPEC. 

We have checked that our 15-T model based on Osten et al.'s model was 
correct by convolving it through the same Chandra response matrix files and 
comparing it with their synthesized model provided by Dr. Osten. 



the current calibration k nowledg e; iii) it was inadequate to 
multiply the EMs of Oste n et all s model by a constant value 
of 1.39 to account for the X-ray luminosity difference of 
CrB between the Chandra and XMM-Newton observations. 
It is possible that the variation is due to some part of the 
DEM only (e.g., at high temperatures). However, the long- 
wavelength continuum level is essentially determined by the 
high-temperature components (which describe the short and 
mid-wavelength ranges accurately), and, therefore we find 
this possibility less probable. Instead we propose that the 
lack of spectral coverage in the range 25-40 A by Cha n- 
dra HETGS sUghtly biased the EMD of Oste n et alJ (120031) to 
overestimate the long-wavelength flux. We note that, although 
Osten et al. (2003) used EUVE line fluxes to determine their 
DEM, these authors did not make use of the continuum level 
in the EUVE wavelength range (which is about 5 times lower 
than the continuum level at 30 A). 

Our XMM-Newton model (whose EMs were divided by 
1.39) sho ws gene rally excellent agreement with the Chandra 
data (Fig. I6al6c> . However, the flux in the Ar XVII triplet 
near 4 A was underestimated; it suggests that Ar abundances 
obtained from simultaneo us EPIC-i-RG S fits might be under- 
estimated generally. The Ost en et all model, however, pre- 
dicts L-shell fluxes of Ar and S lines at long wavelengths 
which are consistent with the RGS data (Fig. ISct . despite 
the long-standing problem of the lack of atomic data of L- 
shell lines of low- Z elements or their i naccuracy in current 
atomic databases ( A udard et a Ll l2001ht iRaassen et aU^XKY^L 
iLepson et al. 2003). 

With a few exceptions, we find, therefore, excellent agree- 
ment among different models derived with different tech- 
niques from different data sets, but which show very simi- 
lar abundance ratios but different absolute abundances. This 
exercise casts some doubt on the preference of one method 
over the other to determine coronal abundances, and it demon- 
strates that the determination of true absolute abundances re- 
mains a chal lenge.^ Other stud ies support this interpretation. 
For example. lAudard et al.'('2004) studied the bright FK Com- 
type giant YY Men using different methods and found that, 
for a given data set (either Chandra or XMM-Newton), slight 
differences in absolute abundances were derived. However, 
they found excellent agre ement for the abundan ce ratios rel- 
ative to Fe. In parallel, Telleschi et alJ (12005^ studied the 
XMM-Newton data of a sample of solar analogs with different 
methods as well. They also observed the robustness of abun- 
dance ratios and demonstrated the presence of a transition 
from a solar-like FIP effect in inactive stars t o an in verse FIP 
effect in active stars. iGarcia- Alvarez et'al] (l2005t) used one 
method but different sets of emission lines to derive the coro- 
nal abundances in AB Dor (and V471 Tau) and found similar 
abundances. They further compared their results with those in 
earlier studies and found good agreement in the general trend, 
albeit with some significant differences. Thus, in brief, de- 
spite the use of various methodologies, abundance patterns as 
a function of the FIP appear robust and look dependent on the 
magnetic activity level. However, it remains challenging i) 

^ Although not shown in this paper, we note that setting the absolute Fe 
abundance to solar photospheric and using the same abundance ratios as pro- 
vided by Osten et al. (2003), it is impossible to obtain a reasonably good fit 
of the Chandra data by simply multiplying the EMs by a constant factor If 
the model is adjusted to match the continuum level, the Fe line fluxes are 
indeed systematically overestimated. It again reinforces the result that Fe is 
underabundant in CrB's corona compared to the solar photospheric Fe 
abundance. 
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to constrain absolute abundances, and ii) to compare coronal 
abundances to their photospheric counterparts. 

5.3. Electron Density 

High-resolution X-ray spectra observed with XMM-Newton 
and Chandra allow us to determine electron densities from 
line ratios of He-like ions. The ratio R of the flux in the 
forbidden (/) line to the flux in the intercombination (/) line 
is indeed density-sensitive. In particular, the O VII triplet is 
most useful in measuring coronal densities of the order of 
lO"* cm"^, and is relatively blend-free (in contrast to Ne IX), 
and well-resolved (in contrast to, e.g., Mg XI and Si XIII for 
the RGS). We have obtained line fluxes for the O VII triplet 
using delta functions convolved through the instrumental re- 
sponse. We list them in Table |3 together with 68% confi- 
dence ranges (Ax^ = 1). The R ratio is equal to 2.81 ±0.11, 
corresponding t o an upper limit t o the electron density of 
4 X lO'*' cm (Porauet et al. 2001), assuming a photoexci- 
tation contribution by the radiation field with T^ff ^ 6,000 K. 

Osten et al. ( 2003 ) and Testa et al. ( 2004) reported for O VII 
an electron density of ~ 2 x 1 0'" cm"-', i.e., consistent with 
ours. These authors and ISanz -Forcada et al. (2003b) quote 
(higher) electron densities from other He-like triplets and Fe 
lines (in the EUV range) which we cannot confirm due to the 
poorer spectral resolution of the RGS (or the lack of wave- 
length coverage). Unfortunately, the N VI triplet is too faint 
to measure any line flux. 

5.4. Opacity 

Our CIE models assume an optically thin plasma, so it is 
necessary to verify that opacity effects play a negligible role. 
We examined the ratio of the Fe XVII spectral lines at 15.01 A 
and at 15.26 A, since this ratio is sensitive to opacity effects 
(the 15.01 A line has a large oscillator strength). The mea- 
sured Fe XVII fluxes are shown in Table0] 

Until recently, theoretical values for the ratio used to 
range from 3.0 to 4.7 (Bhatia & Doschekl \l99% . signifi- 
cantly different from laboratory measuremen ts which sug- 
gest a range from 2.8 to 3.2 (iBrown et al.il998t | La ming et al.| 
12000: Brown et al. 2001). However, recent theoretical work 
has brought the theoretic al ratio clo ser to the laboratory ratio 
JChen & Pradhaa.2002: Chen et al.l2003 - .Fournier & HansenI 
I2005h . We obtain a flux ratio in CrB of 3.02 ±0.15 
i.e., c onsistent with the the ratios derived by lOsten et aP 
(l200l (3.0 1 ±0.26 with HEG and 2.94 ±0.15 with MEG). 
iNess et al.1 12003) derived a slightly lower ratio (2.56 ±0.13), 
probably because of a different treatment of the continuum. 
Nevertheless, they also concluded that this value was close to 
the ratio obtained in laboratory measurements. Thus, the line 
ratio indicates that there are no significant opacity effects in 
the coronae of CrB. 

6. SUMMARY AND CONCLUSIONS 

We have presented a study of the XMM-Newton observa- 
tions of the binary CrB. No large flare was detected dur- 
ing the observations. We focused our analysis on the de- 
termination of the coronal abundances by means of a multi- 
T approach. Although such models do not give an accu- 
rate representation of the coronal DEM, they suffice to obtain 
good measurements of c oronal abund ances and abundance 
ratios in particular (e.g., lAudard e t al."2004: Telleschi et alj 
l2005h . We have, furt hermore, compared our resulting abun- 
dances with those of lOsten et alJ J2003h who published the 



Chandra data of cr^ CrB. We find very good agreement be- 
tween their abundance ratios (relative to the Fe abundance) 
and ours, but find a systematic difference in the absolute 
abundances. Upon checking whether their model could fit 
our XMM-Newton data, we have obtained an excellent agree- 
ment (after correcting the published DEM, see Section I5T2I 1. 
At wavelengths longer than 25 A, their model overestimates 
the continuum level. We argue that their model, based on 
HETGS data below this wavelength, did not constrain their 
DEM well enough to account for the longer wavelength range. 
Our XMM-Newton best-fit model also reproduces the Chan- 
dra HETGS spectra remarkably well. Thus, abundance ra- 
tios are found to be robust, whereas the discrepancy in abso- 
lute abundances found between our model and Osten et al.'s 
model probably stems from the different methodologies. Sim- 
ilar re sults hav e been found elsew here (^Audard et al. 20041 
Telles chi et alJ l2005: Gai'cia-Alvar ez et alJl2005l) . Neverthe- 
less the presence of abundance patterns as a function of the 
FIP seems reliable. 

The abundance ratios relative to the coronal Fe abundance 
showed a correlation with the FIP, but such a correlation was 
neither similar to the inverse FTP eff ect observed in very ac- 
tive binaries (e.g. lAudard et alJ2003^ nor to the solar-like FIP 
effect observed in less active stars (e.g.pelleschi et al. 2005). 
The abundance ratios of low-FIP elements (< 10 eV) with re- 
spect to Fe are similar to the solar photospheric abundance ra- 
tios; for intermediate-FIP elements, they are however consid- 
erably lower but increase with increasing FIP. Such an abun- 
dance patter n has already been observed in other magnetically 
active stars (Gudeletal. 2001a; Sanz-Forcada et al. 2003al 
Huenemoerder et al. 2003). However, Garcia- Alvarez et all 
(2005) found the abundances of the very low-FIP elements 
(Al, Ca, Na) in AB Dor to be similar to, or lower than, those 
of low FIP (Mg, Si, Fe), in contrast to Sanz-Forcada et al.'s re- 
sults. It is possible that these patterns are real, but we empha- 
size that coronal abundances are often compared to the solar 
photospheric composition due to the lack of available or reli- 
able measurements in magnetically active, fast-rotating stars. 
In the case of a~ CrB, the Fe and Ca phot ospheric abundances 
of the CrB binary were available ( Strassmeier & Ric^ 
I2OO3) and suggested that the low Fe abundance in CrB's 
coronae was in fact consistent with the photospheric value. 

Finally, we obtained an upper limit of the coronal density 
based on the O VII triplet (wg < 4 x 10'" c m~^) which is con- 
sistent with other measurements in a CrB (Osten et al. 20031 
^staetal. 2004). Osten et al. (2003), Sanz-Forcada et aQ 
(2003b), and Testa et al. (2004) found higher densities (Wg ^ 
10'^ cm"^) from spectral lines formed at higher temperatures 
than O VII. However, due to the lower spectral resolution of 
the RGS (compared to Chandra HETGS), we could not con- 
firm such coronal densities. 

In conclusion, this study has confirmed that the abundance 
ratio pattern in CrB's coronae does not follow that of a 
simple inverse FIP effect or of a solar-like FIP effect, but ap- 
pears to be more complex. However, the lack of data on pho- 
tospheric abundances in CrB to compare with still casts 
doubt on whether such a pattern is real or mirrors the stellar 
composition. Nevertheless, the underabundance of Fe in <t^ 
CrB's coronae (compared to the solar photospheric Fe abun- 
dance) appears robust and is consistent with its abundance 
found in the stellar photosphere. 
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TABLE 1 
Observation Log. 







Observation Start 


Observation End Total Exposure 


Filtered Expostire 


Instrunieiil 


JVloclc 


FUter (UT) 


(UT) 


(ks) 


(KS) 






Revolution 3 16, Observation ID 0111470201 ("201") 






RGSl . . . . 


Spectroscopy 


None 2001 Aug 31 07:38:02 


2001 Aug 31 13:02:20 


19.5 


18.9 


RGS2 .... 


Spectroscopy 


None 2001 Aug 31 07:38:02 


2001 Aug 31 13:05:11 


19.6 


18.5 


MOSl''. .. 


TIMING 


Thick 2001 Aug 31 07:44:30 


2001 Aug 31 11:32:36 


13.7 




MOS2. . . . 


IM.AGING 


Medium 2001 Au2 31 07:44:33 


2001 Aus 31 1 1:34:12 


13.8 


9.8 


PN 


TIMING 


Medium 2001 Aug 31 08:13:58 


2001 Aug 31 13:08:00 


17.6 








Revolution 404, Observation ID 01 1 1 47030 1 ("301") 






RGSl . . . . 


Spectroscopy 


None 2002 Feb 22 23:23:29 


2002 Feb 23 05:47:04 


23.0 


22.5 


RGS2 . . . . 


Spectroscopy 


None 2002 Feb 22 23 :23 :29 


2002 Feb 23 05:47:05 


23.0 


21.8 


MOSl''... 


TIMING 


Thick 2002 Feb 22 23:29:57 


2002 Feb 23 05:27:24 


21.4 




MOS2. . . . 


IMAGING 


Medium 2002 Feb 22 23:29:59 


2002 Feb 23 05:27:53 


21.5 


6.7 


PN'' 


TIMING 


Medium 2002 Feb 22 23:59:25 


2002 Feb 23 05:46:27 


20.8 




Note. — 


Data from Revolution 315 are not included for analysis due to short exposixre times of approximately 7ks and significant baclcground 


contamination from solar flare activity. 









^Maximum effective integration time used for data analysis. 
''The EPIC MOS 1 and pn data were not used in this paper. 
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TABLE 2 
Wavelength Ranges 
Excluded from the 
Fitting Procedure. 



Instrument X Range 

(A) 



RGS 


< 8.3 




RGS 


15.80-16. 


,20 


RGS 


18.60-18. 


,80 


RGS 


19.50-19. 


,90 


RGS 


20.00-20. 


,30 


RGS 


20.60-20, 


,80 


RGS 


21.15-21. 


,30 


RGS 


21.40-21. 


,50 


RGS 


23.40-23. 


,70 


RGS 


24.10-24. 


,60 


RGS 


24.90-25, 


,25 


RGS 


25.45-25, 


,80 


RGS 


26.25-26, 


,50 


RGS 


27.00-27, 


,20 


RGS 


27.30-27, 


,70 


RGS 


27.80-28, 


,10 


RGS 


28.30-28, 


,50 


RGS 


29.40-29, 


,70 


RGS 


30.20-30, 


,60 


RGS 


30.90-31. 


,15 


RGS 


31.50-31. 


,70 


RGS 


32.10-33. 


,55 


RGS 


33.85-34 


,00 


RGS 


34.50-34. 


,65 


RGS 


35.60-35. 


,80 


RGS 


36.30-36. 


,85 


RGS 


37.50-37. 


,80 


M0S2 . . . 


< 1.65 




M0S2 . . . 


> 18 
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TABLE 3 

Best-Fit Parameters for RGS+M0S2 Togeth er with 
Abundances from IOsten et al.I 120031) . 



Parameters 201 

[C/Fe] -O-^^Tlo 

[N/Fe] ""-"^^U 

[O/Fe] -'^■^^m 

[Ne/Fe] +"■^1^03 

[Mg/Fe] +0-04!^ o| 

[Si/Fe] -0.06:jj 

[S/Fe] -O-^S^nig 

[Ar/Fe] +0.04;[|| 

[Ca/Fe] -0.12;!i;30 

[Fe/H] -0-17+!!m 

(MK) 3.83!|! jO 

7-2 (MK) 7.74!|!U« 

7-3 (MK) 14.7:j! 5^ 

r4(MK) ^^-4^' 13 

logEMi (cm-3) 52.7™:]0 

logEM, (cm-3) 53.l!| ]p 

logEMs (cm-3) ^^-^^ lo 

logEM4 (cm-3) 52.7!«:]0 

1-34 



301 



Osten et al. (2003) 



-0.22 



+0.15 

+o.08!«;»^ 

-0.09!«:i 

-0 49+01'' 
-Q.32 
-0 35+''-^2 

+0 16+"-2^ 
-Q.53 
-0 13+"''2 

T ^4+0. 18 
^- -0.19 
7 4K+0.08 

1 4 4+0.50 

31 6+5 "0 

■'^■■'-o.io 
53 j+0.10 

52 9+S:l8 

0+8:18 

1.48 



-0.07™ 
-0.16™ 
+0.16™ 
-O.OO^S 
-0.10™ 
-0.28+'|' 
+0.34+^ 



-0 34+O 02 

"■-^ -0.02 
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TABLE 4 

Photon Fluxes for O vii 

TRiPLiii' AND Fli xvii Lines. 



Wavelength 


Flux 


(A) 


(10^ s-') 


2L60 


3.65 ± 0.28 


2L80 


0.80 ± 0.20 


22.10 


2.25 ± 0.25 


15.01 


13.9 ± 0.30 


15.26 


4.60 ± 0.21 
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0.36 



0.76 



Phase 
0.86 0.96 



1.06 





1.19 



52150.3 52150.4 52150.5 52150.7 52152.3 52152.4 52152.5 52152.6 52327.9 



52328.1 



52328.2 52328.3 



Time (MJD) 



Fig. 1 . — EPIC MOS2 lightciu-ves with a bin size of 500 s for the three observation dates (error bars are shown). The selected good time intervals are shown by 
horizontal lines above the lightcurves. The excluded time ranges are during high background solar flare activities. The net count rates are reasonably constant. 
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14 16 
Wavelength (A) 



18 



Fig. 2a.— (a) XMM-Newton MOS2, RGSl, and RGS2 spectra with the EPIC MOS2+RGS (301; see Tab.0 best-fit model (superposed line) for 1-20 A. 
Error bars are not shown for clarity. The model is only shown in the wavelength ranges used for the spectral fits (see Tab.l2l. The ranges for the RGS are also 
shown at the top of each panel as broken lines, {h) Same as (a), but for 20—40 A. 
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Fig. 2b.— 
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First Ionization Potential (eV) 



20 



Fig. 3. — Abundance ratios relative to Fe and relative to solar abundances IGrevesse & Sauvall998ft for the 201 and 301 observations as a function of the FIR 
For clarity, we shifted the 301 data points by +0.2 eV. Statistical errors only are shown; systematic errors probably range up to 0. 1 —0.2 dex. For Ca (201) and Ar 
(301), upper limits are shown with an an'ow down to the bottom of the figure. 
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[A/Fe], 301 
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-0.5 0.0 

[A/H], 201 



1.0 



Fig. 4. — Comparisons between abundances obtained in this paper and with those by Osten et al. 1 2003). The bottom right panel (d) shows absolute abundances 
instead of abundance ratios to Fe (a,b,c). Abundance ratios are robust and similar, whereas our absolute (relative to H) abundances are systematically higher 
than those derived by Osten et al. ( 2003). Some abundances discussed in the text are labeled for clarity. The dotted line represents a 1:1 correlation, whereas the 
dashed line in panel d represents a 1.5:1 correlation. 
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